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ABSTRACT 

We analyze the Far-UV/UV spectra of four central stars of planetary nebulae 
with strong wind features — NGC 2371, Abell 78, IC 4776 and NGC 1535, 
and derive their photospheric and wind parameters by modeling high-resolution 
FUSE (Far- Ultraviolet Spectroscopic Explorer) data in the Far-UV and HST- 
STIS and lUE data in the UV with spherical non-LTE line-blanketed model 
atmospheres. Abell 78 is a hydrogen-deficient transitional [WR]-PG 1159 object, 
and we find NGC 2371 to be in the same stage, both migrating from the constant- 
luminosity phase to the white dwarf cooling sequence with T^ff ~ 120 kK, M ~ 
5xl0~^ MqYt-\ NGC 1535 is a "hydrogen-rich" 0(H) CSPN, and the exact 
nature of IC 4776 is ambiguous, although it appears to be helium burning. Both 
objects lie on the constant-luminosity branch of post-AGB evolution and have 
Teff ~ 65 kK, M ^ lxlO~^. Thus, both the H-rich and H-deficient channels of 
PN evolution are represented in our sample. We also investigate the effects of 
including higher ionization stages of iron (up to Fe X) in the model atmosphere 
calculations of these hot objects (usually neglected in previous analyses), and 
find iron to be a useful diagnostic of the stellar parameters in some cases. The 
Far-UV spectra of all four objects show evidence of hot (T ~ 300 K) molecular 
hydrogen in their circumstellar environments. 



Subject headings: stars: atmospheres — stars: Wolf-Rayet — stars: individual 
(NGC 2371, Abell 78, NGC 1535, IC 4776) — planetary nebulae: general — UV: 
spectroscopy 



"'^Based on observations made with the NASA-CNES-CSA Far Uhraviolet Spectroscopic Explorer and 
data from the MAST archive. FUSE is operated for NASA by the Johns Hopkins University under NASA 
contract NAS5-32985. 
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1. INTRODUCTION 

After leaving the asymptotic giant branch (AGB), low and intermediate- mass stars 
become central stars of planetary nebulae (CSPN). The spectral types of the central stars 
are numerous and diverse {e.g., O, Of, sdO, [WR], [WR]-PG 1159, as well as white dwarf- 
types), but all of these can be generally divided into two groups based on surface abundances: 
hydrogen-rich (which show obvious hydrogen lines in their spectra) and hydrogen-deficient 
(which don't). H-deficient stars are in a post-hehum flash, helium-burning phase, while the 
H-rich stars may be either hydrogen or helium burning. About 10-20 % of CSPN are H- 
deficient (De Marco & Soker 2002; Koesterke & Werner 1998 and references therein). The two 
groups are thought to represent different "channels" of CSPN evolution, terminating with 
either H-rich DA type or He-rich DO type white dwarfs {e.g., Napiwotzki 1999). During the 
CSPN phase the star sheds mass in the form of a stellar wind as it evolves toward the white 
dwarf cooling sequence. It is generally believed that the winds of these objects are driven by 
radiation pressure, as with population-I 0-stars and WR stars. If the driving of the wind is 
mainly dictated by the stellar luminosity, the wind parameters can, in principle, be used to 
determine the distances to CSPN (Kudritzki et al. 1999; Kudritzki & Puis 2000; Tinkler & 
Lamers 2002). 

The spectra of the majority of the H-deficient class have strong wind signatures, show 
no photospheric hnes and are very similar to those of WC stars — the carbon-rich Wolf- 
Rayets which are evolved massive stars. The low-mass central stars are termed "[WC]" stars 
to denote their association with planetary nebulae. "PC 1159"- type stars are extremely 
hot, H-deficient white dwarf-type objects, thought to represent an entry point onto the 
white dwarf cooling sequence. They show mainly absorption lines in the optical spectra but 
possess a few wind lines in the UV. Because PG 1159 objects have similar abundances to 
[WC] stars and about half of the PG 1159 stars are CSPN (Werner 2001), they are thought 
to be descendants of the [WCjs. Linking these two classes is a rare group of objects named 
"[WC]-PG 1159" stars, whose spectra show both emission lines and some absorption lines. 
Only three definite cases of such stars are known so far: Abell 78, Abell 30 and Longmore 
4 (Koesterke 2001). Parthasarathy et al. (1998) have posited that all "weak emission line 
stars" ([WELS]) are actually [WC]-PG 1159, based on optical spectra, but Werner (2001) 
cautioned against this statement until better spectra of the sample is available. The H-rich 
CSPN sample of Mcndez et al. (1988) are similar in that their spectra also feature absorption 
lines as well as UV wind emission features. 

[WC], [WC]-PG 1159, and PG 1159 stars all reside in the same area of the H-R di- 
agram (HRD). In terms of evolution, this region corresponds to the bend where the star 
leaves the constant-luminosity post- AGB phase and is transitioning to the white dwarf cool- 
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ing sequence. [WC] and PG 1159 stars lie close to each other in the H-R diagram, yet 
have very different spectra (in the optical, a [WC] spectrum is littered with strong emission 
lines, while that of a PG 1159 presents an absorption line spectrum). This suggests a dra- 
matic drop in mass-loss rates over a short time of evolution (Koesterke & Werner 1998). 
The commonly accepted evolutionary sequence for H-deficient GSPN is: [WC] — i^[WC]- 
PG 1159— >PG 1159— >D0 WD (Hamann 1996). This sequence spans the time where the 
stellar winds are in the process of "turning-off " . Therefore, understanding the objects in 
this part of the H-R diagram is not only a prerequisite for understanding the evolution of 
central stars, but will also give insight into the driving of stellar winds in general. 

In the classification of [WR]-type objects, wind parameters are more fundamental dis- 
criminators than stellar temperature (Crowther 1999; Acker & Neiner 2003). As the central 
star's winds die down, the wind features in the FUV/UV are the last to fade (Koesterke et al. 
1998). Therefore, to estimate accurate parameters of the extended atmospheres this phase 
is best studied in the FUV/UV. In this paper, we present a FUV/UV-based spectral anal- 
ysis of four CSPN which show strong wind signatures - NGC 2371, Abell 78, IC 4776, and 
NGC 1535. We model the FUV/UV spectra of these stars using stellar atmospheres codes to 
determine the wind parameters such as T^ff, M, and v^o and discuss their evolutionary im- 
plications. Abell 78 has been classified as a [WO]-PG 1159 transition object (Mendez 1991; 
Crowther et al. 1998), and we argue that NGC 2371 also belongs to this class. NGC 1535 
is an 0(H) star belonging to the H-rich sample of Mendez et al. (1988). The FUV spectra 
of IC 4776 is similar to that of NGC 1535, but we find evidence that this object may be 
a He-burner. Comparing with evolutionary calculations, our derived parameters place the 
latter two objects along the constant-luminosity path of the HRD, while the former two lie 
on the bend marking the transition from that phase to the white dwarf cooling sequence. 

This paper is arranged as follows. The observations and data reduction are described 
in § 2. A comparison of the spectra of the objects is presented in § 3. Our models and 
parameter determinations are described in § 4. The implications of our results are discussed 
in § 5 and our conclusions in § 6. 

2. OBSERVATIONS AND REDUCTION 

The coordinates, radial velocities (vr) and nebular diameters (9) of our four sample 
objects are listed in Table 1 (Tylenda et al. 2003 has measured two angular dimensions for 
three of the sample PN, corresponding to the semi-major and minor diameters, which we list 
as "^B X 9a"). The data sets utilized in this paper are summarized in Table 2. NGC 2371 
and IC 4776 were observed as part of FUSE cycle 1 program P133 (Bianchi). All other 
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FUV/UV data were retrieved from the MAST archive, including FUSE and STIS archive 
data of NGC 1535, a Berkeley Extreme and Far-UV Spectrometer (BEFS) archive spectrum 
of Abell 78, and lUE data of most objects. The reduction of the FUSE data is described in 
§2.1. 

2.1. FUV Data 

FUSE covers the wavelength range of 905-1187 A at a spectral resolution of < 30,000. 

It is described by Moos et al. (2000) and its on-orbit performance is discussed by Sahnow 
et al. (2000). FUSE collects light concurrently in four different channels (LiFl, LiF2, SiCl, 
and SiC2). Each channel is recorded by two detectors, each divided into two segments (A & 
B) covering different subsets of the above range with some overlap. 

The FUSE spectra were taken using the LWRS (30"x30") aperture. These data, taken in 
"time-tag" mode, have been calibrated using the most recent FUSE data reduction pipeline, 
efficiency curves and wavelength solutions (CALFUSE v2.2). We combined the data from 
different segments, weighting by detector sensitivity, and rebinned them to a uniform dis- 
persion of 0.05 A (which is probably close to the actual resolution since the data were taken 
in the early part of the mission). Bad areas of the detectors, and those regimes affected by 
an instrumental artifact known as "the worm" (FUSE Data Handbook vl.l), were excluded. 
The flux calibration accuracy of FUSE is < 10 % (Sahnow et al. 2000). 

The BEFS spectra covers a wavelength range of 912-1218 A with a resolution of ~ 
0.33 A. The cahbration and performance are described by Hurwitz (1998). 

2.2. UV Data 

For most stars, we make use of lUE data, using high-resolution spectrum if available (a 
high-resolution spectra for NGC 2371 exists, but is of poor quality, so we present the low- 
resolution data instead). For some stars, there was some disagreement (< 30 %) between 
the FUV (both FUSE and BEFS) and lUE continuum levels in the regions of overlap (flux 
cahbrations for these instruments are typically accurate to ~ 10 %). Because there are no 
nearby objects that could contribute extra flux in the FUSE aperture, we suspect that the 
background has been over-subtracted in the lUE spectra since nebular continuum emission 
is expected (in the case of NGC 1535, the STIS spectra agree perfectly well with the FUSE 
fluxes, but the lUE fluxes arc a factor of ~ 1.2 lower). Therefore we we have scaled the UV 
flux levels to those of the FUV (scaling factors are listed in Table 2). 
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3. DESCRIPTION OF THE SPECTRA 

Fig. 1 displays the combined FUV/UV spectra of our sample (up to 1700 A, longwards 
of this, the spectra are mostly featureless). Much of the fine structure between 900 — 1200 A 
is caused by sight-line molecular hydrogen (H2), and is discussed in § 4.2. As far as stellar 
features go, the spectra present a sparse landscape. The O VI AA1032-38 feature is most 
prominent in all cases, and the spectra are arranged from top to bottom in order of decreasing 
strength of this feature. IC 4776 and NGC 1535 both display S VI AA933-44 in the FUV, 
however, their UV spectra are different. In the FUV, Abell 78 and NGC 2371 have 3 similar 
lines but N V AA1238-43 and O V A1371 are strong P-Cygni in Abell 78 while absent in 
NGC 2371, neither showing S VI but both displaying a P-Cygni feature around 977 A. 
This feature is commonly assumed to be C III, however the formation of a C III line in the 
highly-ionized winds of these objects is hard to explain, and this will be discussed later (§ 5). 

C IV AA1548-51 appears in all spectra, although it is relatively weak in NGC 1535. N V 
AA1238-43 is obviously present in all spectra save that of NGC 2371. O V A1371 is absent 
in NGC 2371, present as a P-Cygni profile in Abell 78 and NGC 1535, and appears perhaps 
only in absorption in IC 4776. He II A1640 appears to be present in the low-resolution 
lUE spectrum of NGC 2371, but inspection of the one available (noisy) high-resolution lUE 
spectra for this object (SWP53006) shows the He II A1640 to be at least partly nebular in 
origin. The absence of transitions from ions of lower ionization potential {e.g., N IV, O IV, 
or Si IV) indicate the high ionization state of the winds, and high effective temperatures of 
the central stars. 



4. MODELING 

The modehng of these spectra consists of three parts: modeling the central star spec- 
trum, the nebular continuum emission longwards of 1200 A, and the sight-line Hydrogen 
(atomic and molecular) absorption shortwards. The reddening must also be adjusted con- 
currently to match the FUV/UV slope. In all four objects, the reddening is low, and not a 
significant source of uncertainty. Iterations between each are required to achieve a consistent 
solution. However, we discuss each in turn for clarity. 

4.1. The Central Stars 

The optical spectrum of these objects shows mainly absorption lines, and some have 
been modeled using plane-parallel codes {e.g., Abell 78 by Werner et al. 2003, NGC 1535 
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by Mendez et al. 1988). In the FUV/UV, the most prominent features are wind hnes. The 
temperatures and wind parameters of some of these objects have been previously determined 
from plane-parallel analyses in conjunction with wind-line analyses of one or two wind fea- 
tures (such as C IV AA1548-51 or O VI AA1332-38 — Koesterke & Werner 1998). Here we 
derive wind parameters solely from the FUV/UV spectra. 

Intense radiation fields, a (relatively) low wind density, and an extended atmosphere 
invalidate the assumptions of thermodynamic equilibrium and a plane-parallel geometry for 
the winds of our sample CSPN. To model these winds, we use the non-LTE (NLTE) line- 
blankctcd code CMFGEN of Hillicr & Miller (1998. 1999). CMFGEN solves the radiative 
transfer equation in an extended, spherically-symmetric expanding atmosphere. Originally 
developed to model the Wolf-Rayet winds, it has been adapted for objects with weaker 
winds such as 0-stars and CSPN as described in HiUier et al. (2003). The detailed workings 
of the code are explained in the references therein. To summarize, the code solves for the 
NLTE populations in the co-moving frame of reference. The fundamental photospheric/wind 
parameters include T^ff, M. the elemental abundances and the velocity law (including 
Voo). The stellar radius (i?*) is taken to be the inner boundary of the model atmosphere 
(corresponding to a Rosseland optical depth of ~ 20). The temperature at different depths 
is determined by the stellar temperature T^, related to the luminosity and radius by L = 
47ri?*^(7T'*^, whereas the effective temperature (Teff) is similarly defined but at a radius 
corresponding to a Rosseland optical depth of 2/3 (-R2/3)- The luminosity is conserved at 
all depths, so L = An R2/3^crTgff'^. 

The terminal velocity (voo) can be estimated from the blue edge of the P-Cygni ab- 
sorption features, preferably from features that extend further out in the wind (we used 
C IV AA1548-51 primarily). We assume what is essentially a standard velocity law v{r) — 
Voo(l — To/r)^ where Tq is roughly equal to i?*. Once a velocity law is specified, the density 
structure of the wind is then parameterized by the mass-loss rate through the equation of 
continuity: M — An Rjp{r)v{r). 

It has been found that models with the same transformed radius Rf [oc i?*(voo/ M)^/^] 
(Schmutz et al. 1989) and Vqo have the same ionization structure, temperature stratification 
(aside from a scaling by R^,) and spectra (aside from a scaling of the absolute flux by R^^ 
— Schmutz et al. 1989; Hamann et al. 1993). Thus, once the velocity law and abundances 
are set, one parameter may be fixed (say i?*) and parameter space can then be explored 
by varying only the other two parameters {e.g., M and T^ff). Rt can be thought of as an 
optical depth parameter, as the optical depth of the wind scales as oc Rt~^, for opacities 
which are proportional to the square of the density. Scaling the model to the observed flux 
yields R^/D — a distance D must be adopted to determine i?*. 
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For the model ions, CMFGEN utilizes the concept of "superlevels" , whereby levels of 
similar energies arc grouped together and treated as a single level in the rate equations 
(Hillier & Miller 1998). Ions and the number of levels and superlevels included in the model 
calculations, as well references to the atomic data, are given in the Appendix (§ A). 

4.1.1. Clumping 

Radiation driven winds have been shown to be inherently unstable (Owocki et al. 1988, 
1994), which should lead to the formation of clumps in the expanding atmosphere. The 
degree of dumpiness is parametrized by /, the filling factor. One actually can only derive 
Ms={ Mc/y/f) from the models, where Ms and Mc arc the smooth and clumped mass-loss 
rates. For the denser winds of population-1 Wolf-Rayet stars, the clumping factor can be 
constrained by the strength of the electron scattered optical line wings, and clumping factors 
of / ~ 0.1 are typical (a reduction of M to roughly a third of its smooth value). For 0-stars, 
the lower mass loss rates make the electron scattering effects small and difficult to constrain 
(Hillier et al. 2003). The winds of our sample stars are even weaker, and we do not attempt 
to determine the clumping for our sample. Unless otherwise noted, we have adopted f — 1 
and use M to refer to the smooth mass loss rate throughout this paper. This can be seen as 
an upper limit. One expects the clumping in the more tenuous CSPN wind to be less severe 
than in a WR wind, so a conservative estimate of the lower limit of M would be a third of 
this value. 

4.1.2. Gravity 

Because of the FUV/UV range contains mostly wind lines, there are no suitable ab- 
sorption lines to be used as gravity diagnostics. In CMFGEN, gravity enters through the 
scale height h (oc g~^), which connects the spherically extended hydrostatic outer layers to 
the wind. The relation between h and g, defined in Hillier et al. (2003), involves the mean 
ionic mass and mean number of electrons per ion, the local electron temperature and the 
ratio of radiation pressure to the gravity. Mendez et al. (1988) determined logg' = 4.3 ± 0.2 
for NGC 1535 and Werner et al. (2003) found logg = 5.5 for Abell 78 through analysis of 
photospheric optical lines. We initially adopted log g — 4.3 for NGC 1535 and IC 4776, and 
logg = 5.5 for Abell 78 and NGC 2371. For such gravities, the derived wind parameters are 
not that sensitive to h. Spectral wind features can be fit with the models of the same T^ff 
but with gravities differing by over a magnitude. We discuss the issue of gravity more in § 5. 
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4.1.3. Abundances 

Throughout this work, the nomenclature Xi represents the mass fraction of element i, 
"Xq" denotes the solar abundance, and our values for "solar" are taken from Gray (1992). 

Because of the lack of photospheric lines, we have made assumptions about the abun- 
dances of the model atmospheres. We have computed models with two different abun- 
dance patterns: solar abundance for all elements (corresponding to H-rich objects) and 
carbon/oxygen enriched "[WC]-PG 1159"-type abundances, typical of H-deficient CSPN. 
With respect to the latter, it has been shown that abundances through the [WC] subclasses 
are about the same (Crowther et al. 2002, 1998), with the spectroscopic differences mainly 
tied to M and Te// (Crowther et al. 2002). Calculations by Herwig (2001) show that these 
abundances may result when the star experiences a late He-shell flash. The abundances 
of PC 1159 objects are similar (Gorny & Tylenda 2000). Typical measured values are (by 
mass): XHe= 0.33-0.80, Xc= 0.15-0.50, Xo= 0.06-0.17 (Werner 2001). For these rea- 
sons, we have assumed a similar abundance pattern for these elements {Xue/ ^c/ Xo= 
0.54/0.36/0.08), while adopting solar abundances for others {e.g., S, Si, & Fe). The nitrogen 
abundance in such objects typically ranges from undetectable to ~ 2% (by mass, see Werner 
2001 and references therein). We have adopted a nitrogen abundance of 1% for most of our 
models. The lone N V AA1238-43 feature is abundance sensitive, and we do not rely on it 
strongly as a diagnostic as it is also sensitive to T^ff and M. We use these abundances in 
our models of AbeU 78 (a [WR]-PG 1159 star) and of NGC 2371, which has a similar FUV 
spectrum (Fig. 1). 

For the H-rich 0(H) star NGC 1535, we have adopted solar abundances for all elements. 
For IC 4776, which displays a similar FUV/UV spectrum, we have computed models with 
both solar abundances and the above H-deficient "[WC]/PG 1159" abundances. We are able 
to fit its spectrum adequately with models of both abundances (we will discuss this more in 
§4.1.6). 

Table 3 shows our final model elemental abundances. For most, our assumed abundances 
produced adequate fits. Some required adjusting, as described in § 4.1.6 in a case-by-case 
basis. 

4.1.4. Diagnostics 

Optimally, one determines Te// by the ionization balance of He or the CNO elements. 
It is desirable to use diagnostics from many elements to ensure consistency. However, in the 
FUV/UV spectra of these objects only oxygen shows two ionization stages - O V and O VI. 
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We find the O VI feature to be mainly insensitive to M and moderately sensitive to T^ff, 
however, the possible influence of X-rays on this line make it less than desirable temperature 
diagnostic (in the 0-stars regime, O VI is also sensitive to M, but above a sharp threshold, 
and this line is primarily created by X-rays — Bianchi & Garcia 2002). The other elements 
(C, N, S) only show one ionization stage. Because of the lack of multiple ionization stages 
from the same element in the spectra of most of our sample, effective temperatures are 
mainly derived from the appearance or absence of ions of different ionization potentials. For 
example, in our cooler objects (IC 4776 and NGC 1535) the absence of low ionization features 
such as S V(73 eV), P V(65 eV), O IV(77 eV), and N IV(77 eV) places a lower limit on Teff 
(~55-60 kK from our models) for a wide range of mass-loss rates and abundances. Tinkler 
& Lamers (2002) find that between 55-70 kK, H-rich CSPN exhibit a large jump in M by a 
factor of 10-100, which they suspect is due to a bi-stability jump. In our models, it is in this 
temperature range that these previously mentioned ions jump to the next stage, supporting 
the idea of Vink et al. (2001) that a bi-stability jump at such high temperatures may be 
due to CNO elements. The presence of S VI in these object's spectra further constrains the 
temperature to be below a certain point. Above Teff > 90 kK, S VII becomes dominant in 
the outer wind, and a regime is entered in which our "hotter" objects lie (NGC 2371 and 
Abell 78). Around Te//~120 kK, the winds become too ionized for O V, a discriminator 
between our two hot objects. 

Saturated features (like O VI) are mainly insensitive to the mass-loss rate. We typically 
use the observed C IV, O V, and S VI lines to constrain M. Also, He II A1640 (not observed 
in most of our sample objects) is used to place upper limits on this parameter. In some 
cases, we also use the iron spectrum as an additional diagnostic (§ 4.1.5). 

4.1.5. Iron 

The effects of iron (and other line- blanketing elements) are important to consider in the 
modeling of hot stars for multiple reasons. It is known that in the modeling of hot stars, 
the neglect of important line-blanketing elements such as iron can significantly impact the 
derived parameters. Furthermore, iron features may be used to provide additional diagnostics 
of stellar parameters {e.g., Herald et al. 2001). Finally, line-blanketing elements such as iron 
are high sources of opacity in the wind and are thus important in the understanding of how 
it is driven. 

Our sample represents an ideal test-bed for studying the effects of the inclusion of higher 
ionization stages of iron (i.e., Fe VII-Fe X) in the model atmospheres, which have typically 
been neglected in previous studies of similar objects. The temperatures of the objects in 
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our sample span a range where different stages in tfiis sequence dominate. Since Fe IV-Fe X 
liave transitions wliicli occur in tfie FUV and UV, tfie wavelengtfi coverage of our spectral 
data is adequate to study these effects. 

The ionization structures of iron for models in the temperature range of our sample 
(with similar mass-loss rates) are shown in Fig. 2. The figures illustrate what ionization 
stages are dominant in various temperature regimes. For the cooler objects, Fc VI-Fe VII 
dominate. As T^ff increases, Fe VIII becomes dominant throughout the wind (between 
Tgff > 100 kK) and at the highest temperatures (Tg// > 120 kK) Fe IX becomes important 
in the inner winds layers. Useful iron diagnostics include the array of Fe V transitions 
occurring between 1350-1500 A, and the forest of Fe VI-Fe VII lines between 1250 and 
1500 A. Fig. 3 demonstrates how the iron spectrum can be used as a diagnostic of wind 
parameters for NGC 1535. 



4.1.6. Stellar Modeling Results 

Our model fits arc shown in Fig. 4. We now discuss our results for the individual objects, 
along with some notes. Our derived stellar parameters and a global discussion is given in 
§5. 

NGC 1535: The strongest wind features (S VI, O VI, N V, O V, C IV) of this 0(H) star 
are all well fit by our model using solar abundances. As discussed in § 4.1.5, the iron 
spectrum provides an additional constraint for this 0(H) star, as at higher temperatures the 
Fe VII features are too strong (reducing M serves to weaken the non-iron wind features to 
unacceptable levels). 

Our derived parameters (Te// = 66 ± 5 kK, log M = —8.10 ± 0.3 yr^^) agree 
within the uncertainties with those of Koesterke et al. (2004), who model this object using a 
different wind code {Tefj = 70 kK, log M — —7.8 Mq yr"-*^). The temperature is about the 
same as that found by Mendez et al. (1988) in their plane-parallel analysis of optical spectra 
{Teff = 58 ± 5 kK, logg = 4.3 ± 0.2). It is cooler than that derived by Tinkler & Lamers 
(2002) who reported the following parameters for NGC 1535: T^// = 80 kK, = 0.38 R©, 
logL = 3.76 Lq, log^ = 5.05, v^o = 1900 kms-\ and log M = -8.7 ± 0.7 M© yr'^ 
{Rt = 430l28o)- Their mass-loss rate, which is about a factor of 4 lower than ours, was 
derived solely from O V A 1371, which we find is quite sensitive to T^ff in this temperature 
range. 

IC 4776: The nebular continuum for this object strongly contributes to the observed UV 
spectrum at longer wavelengths (§ 4.3). We are able to achieve fits of about the same quality 
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for this object using either our H-rich or H-deficient abundances by adjusting Tg// within 
the range 55-70 kK and log M within —7.85 ±0.30 Mq yr~^, with the H-rich models yielding 
slightly better fits (it is our H-rich model which we present in the figures). O VI 1032-38 
is a bit weak at Tg// ~ 60, but fits well at Tg// = 70 kK. The S VI 933-44 absorption 
is too strong in our presented model. It is better fit with the lower limit of our mass-loss 
rate, but the other wind features are then a bit weak. Our models of either H-rich and 
H-deficient abundances overproduce O V A1371, which is weak at best in the observations. 
Cooler models bring in (unobserved) O IV features and lower mass-loss rates weaken all other 
features. Reducing the oxygen abundance to 1/lOOth solar brings the feature in fine with 
observations, but then the O VI doublet in the FUV range is far too weak. Bouret et al. 
(2003) had difficulties fitting this feature in O dwarf stars, which have mass-loss rates of 
order 10"^ — 10^^ Mq yr^^ and are cooler ( Tg/j ~ 40 — 50 kK). They found a very clumped 
wind (/ = 0.01) alleviated the problem somewhat. 

Our mass-loss rates also agree with the predictions we obtained using the Vink et al. 
(2001) prescription (§ 5), which uses normal (i.e., non-He or C enriched) abundances. The 
low mass-loss rate of these objects and the dearth of diagnostics make the abundances 
difficult to constrain, and from our analysis, we cannot tell if IC 4776 is H-rich or H-deficient. 
Comparison with the evolutionary tracks (also see § 5) indicate this objects is most probably 
a He-burner. 

IC 4776 has been classified as a [WC6] star (Fcibclman 1999) and as a "weak emission 
lined object" (a WELS — Tylenda et al. 1993). The [WC6] classification for IC 4776 was 
based on the equivalent width of C III A2297. This classification is questionable for the 
following reasons. The spectra of a [WC5] or [WC6] object should contain relatively strong 
C III features throughout its spectrum {e.g., Crowther et al. 1998; Hillier & Miller 1999), 
whereas the FUV/UV spectrum of IC4776 is relatively featureless. The FUSE spectrum of 
this object shows high ionization features such as O VI and S VI, and indicating a state of 
higher ionization than one would expect in a [WC6] object. 

Ahell 18: The prominent wind fines of this object (O VI, N V, O V, C IV) are all well- 
reproduced by our model, with the exception of the strong P-Cygni feature at the location 
of C III A977. Koesterke & Werner (1998) also had problems reproducing this feature in 
Abell 78, and speculated that neglected iron lines might sufficiently cool the outer layers 
of the (otherwise hot) wind to allow for the formation of C III. However, we have tested 
models accounting for ionization stages of iron up to Fe X and had no success (the dominant 
ionization stages in the model winds for these objects are Fe VIII and Fe IX). As iron is 
likely depleted in this object (see below), this is probably not the solution (although cooling 
from other species is possible). 
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Using a solar value for the iron abundance ([Fe]=7.67), we could not simultaneously 
fit the iron spectrum and our other diagnostic lines. For any appreciable mass-loss rate, 
relatively strong (unobserved) iron P-Cygni features are seen in models with temperatures 
spanning ~70 — 130 kK. For lower temperatures ( T^ff < 95 kK), the set of Fc VII lines shown 
in Fig. 3 are seen. The hotter models display strong Fe VIII and/or Fe IX features. This 
supports the findings of Werner et al. (2003), who found Abell 78 to have an iron abundance 
of at most ~0.03 solar. Using this iron abundance, the model spectrum agrees with the 
observations. Miksa et al. (2002) have also found iron deficiency in a large sample of PG 1159 
stars, the supposed descendants of these transitional objects. Iron deficiencies in these 
objects may result when material in the He-intershell is exposed to s-process nucleosynthesis 
during a thermally pulsating AGB or post-AGB phase (Lugaro et al. 2003; Herwig et al. 
2003). 

Koesterke (2001) determined = 115 kK, Rt = 45.2 and Voo = 3750 kms~^ for 
this object, with mass fractions for C/O/N of 40/15/2 %, respectively. Our parameters are 
essentially the same (36/8/1), except we derive Voo = 3200 ± 50 from the fitting the C IV 
AA1248-51. Werner et al. (2003) derived Tg// = 110 kK and log^ = 5.5, also in hue with 
our determinations. 

Abell 78 is the prototype " [WC]-PG 1159" star (Hamann 1996; Koesterke & Werner 
1998) and was originally placed in a transition phase between the [WR] and PG 1159 [i.e., 
0(C)] phases by Mcndcz (1991). Crowthcr et al. (1998) suggested a [W01]-PG 1159 classi- 
fication based its high O VI A3818/C IV A5808 ratio. This object also falls into the category 
termed "O Vl-type CSPN", distinguished by Smith & Aller (1969) from the [WC] class for 
having O VI AA3811-34 as one of their most prominent optical features. Our model param- 
eters (derived solely on the FUV and UV spectral region) adequately reproduces the optical 
wind features of Abell 78. 

Abell 78 is also of interest because it belongs to a small group of PN which have H-poor, 
dusty ejecta. These nebulae consist of a H-poor shell surrounded by a an outer H-rich region, 
indicating a 2nd ejection event of H-deficient material during the post-AGB phase (Medina 
& Pena 2000). 

NGC 2371 : As can be seen in Fig. 1, the FUV flux of this object is almost identical to that 
of Abell 78, except the wind lines are more conspicuous and broader, indicating a higher 
terminal velocity and mass-loss rate. The UV spectrum lacks O V, which, in our models, 
disappears for T^ff >115 — 120 kK. There is a weak signature of N V 1238-43. Because of the 
low resolution of our data, it is unclear if its origin is nebular or stellar. Assuming the latter, 
this feature weakens enough to fit the observations in our synthetic spectra for models with 
Tgff > 125 kK (it should be kept in mind we have assumed a mass fraction of Xn = 0.01 
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for nitrogen). The iron spectrum further constrains the temperature to be Te// ^ 130 kK, 
because the Fe IX features in the FUV are too strong for lower T^ff. However, if one assumes 
that this object is iron- deficient as is Abell 78 (see § 4.1.6), this restriction is removed. As 
with AbeU 78, our models fail to reproduce the strong wind feature at 977 A. 

The low resolution lUE spectrum shows He II A1640, but the one high resolution lUE 
spectrum available shows that at least part of the He II emission is nebular. Our upper limit 
to the mass-loss rate assumes it is all stellar, and our lower limit assumes no stellar He II 
line. The high-resolution spectrum also reveals the C IV emission to be partly nebular. The 
absorption trough appears to be saturated (the high-resolution spectrum is under-exposed, 
and hence extremely noisy, but around the stronger lines there is enough flux to see these 
details) . 

Tylenda et al. (1993) classified NGC 2371 as a [WC3]. However, we find the ionization 
of the wind too high for such a classification. Our model parameters for NGC 2371 are 
close to those of Abell 78 (PG 1159-[WOl]), with the former star being a bit hotter with a 
denser wind. Their synthetic optical spectra are qualitatively similar, and so we suggest a 
[WO]-PG 1159 classification for NGC 2371, but at a higher state of ionization marked by 
the absence of O V and the presence of O VI. This classification needs to be confirmed with 
optical spectra, however. 

Kaler et al. (1993) noted a pecuharity with the optical O VI emission features of this 
star: they consisted of two components, a broad, blended feature (with Voo— 3400 kms~^) 
and two resolved narrower features. It is the only known O VI PNN to display both the 
broad and narrow O VI features. They believed these narrow features were associated with 
the wind rather than the nebula, because their long-slit CCD spectrum showed no extension 
of the narrow O VI features beyond the profile of the CS, confining the O VI formation zone 
to have a radius of ~2000 AU. It appears the O VI AA1032-38 for NGC 2371 also has this 
two component structure, as shown in Fig. 8. Close inspection reveals a narrow emission 
line superimposed on top of the broad P-Cygni profile, which could correspond to the blue 
component of the O VI doublet (at ~ 20 km s"^ in the rest wavelength of the star). The red 
feature may also be present, but H2 absorption prevents a firm conclusion. These narrow 
O VI components may be evidence of shocked material in the nebula. Herald & Bianchi 
(2004) found similar (seemingly nebular) O VI features in the LMC CSPN SMP LMC 62. 

In their study of O VI PNN, StangheUini et al. (1995) have found a correlation between 
the strength of O VI AA381 1-3834 and stellar luminosity. According to the models of our two 
transition stars, this holds true. Throughout our sample, the O VI 1032-38 feature grows 
stronger with increasing effective temperature and luminosity. 
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4.2. Molecular and Atomic Hydrogen 

Absorption due to atomic and molecular Hydrogen (H2) along the sight-line complicates 
the spectra of these objects in the FUSE range where numerous H2 transitions from the 
Lyman (i?^S+-X^S^) and Werner (C^n^-X^S+) sequences lie. Toward a CSPN, this 
sight-line Hydrogen typically consists of interstellar and circumstellar components. Material 
comprising the circumstellar H I and H2 presumably was ejected from the star earlier in its 
history (during the AGB phase), and is thus important from an evolutionary perspective. 
If hot {e.g., T > 300K), a small column density [e.g., log(A'") ~ 16 cm~^] of H2 can lead 
to a complex absorption spectrum, obscuring the underlying stellar spectrum (see Herald 
& Bianchi 2004, Fig. 5). It is therefore necessary to model the H I and H2 characteristics 
to discern the underlying stellar spectrum. Our parameter determinations for sight-line H I 
and H2 are listed in Table 4. 

H2 absorption effects were applied to the model spectrum in the following manner. For 
a given column density (A^) and gas temperature (T), the absorption profile of each line 
is calculated by multiplying the line core optical depth (tq) by the Voigt profile H{a,x) 
(normalized to unity) where x is the frequency in Doppler units and a is the ratio of the 
line damping constant to the Doppler width (the "b" parameter) . The observed flux is then 
Fobs — exp [—ToH{a, x)] x Fintrinsic- To fit the H2 spectrum of a given object, we first as- 
sume the presence of an interstellar component with T = 80 K (corresponding to the mean 
temperature of the ISM — Hughes et al. 1971) and Vturb— lOkms"^. The column density 
of this interstellar component is estimated by fitting the strongest transitions. If absorption 
features due to higher-energy H2 transitions are observed, a second, hotter (circumstellar) 
component is added (see example in Fig. 5). The temperature of the circumstellar compo- 
nent can be determined by the absence/presence of absorption features from transitions of 
different J states, and the column density by fitting these features. Iteration between fit- 
ting the interstellar and circumstellar components are performed, as both contribute to the 
lower-energy features. We note that our terminology of "circumstellar" and "interstellar" 
components is a simplification, and basically indicate a "cool" component (assumed inter- 
stellar) and "hot" component (assumed circumstellar). However, the column density derived 
for the cooler "interstellar" component may also include circumstellar H2. 

In the case of NGC 1535, the H2 has been previously modeled by Bowers et al. (1995) 
using ~ 3 A resolution Hopkins Ultraviolet Telescope. They found the data to be well 
fit employing either a one-component model of T = 300 K or a two component model 
with T = 144 and 500 K. Their upper limit for the circumstellar H2 column density was 
logA^(H2) = 18.4 cm~^. This is consistent with our findings, in which we have fit the H2 
absorption (see § 4.2) using two components, a cool component [T = 80 K, logA^(H2) = 
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18.7 and a hot component [T = 400 K, log A^(H2) = 16.7 cm-^], thanks to the higher 

resolution of the FUSE spectra. 

H I column densities were determined in a similar fashion by fitting the profiles of the 
Lya and/or Ly/3 features (our primary diagnostic feature is indicated in Table 4) . For 
Abell 78 and NGC 1535, the available high- resolution data of the Lya hne enable the sight- 
line H I to be constrained tightly, as illustrated in Fig. 6. For IC 4776 and NGC 2371, the 
low-resolution lUE spectra are inadequate for this purpose. For IC 4776, we use the blue 
side of the Ljj3 profile, but for NCC 2371, the P-Cygni absorption from the O VI doublet 
extends far enough blucward that Ljf3 is obscured. Wc therefore simply assumed an H I 
column density of 21.0 for this object for fitting purposes. 

4.3. Nebular Continuum 

Nebular characteristics, taken from the literature, are shown in Table 5. Also presented 
are nebular radii (r„efe) and dynamic (kinematic) ages Tdyn = rneb/vexp, calculated using the 
angular sizes from Table 1 and the distances from Table 7. Since some nebulae have angular 
diameter measurements in two dimensions, we used the average, listed as Oadopt- 

To determine if the nebula significantly contributes to the UV fiux, nebular continuum 
emission models were computed accounting for two-photon, H and He recombination, and 
free- free emission processes. The continuum model parameters are the electron density (rZe), 
the electron temperature (Tg), the Helium to Hydrogen ratio (He/H) and the doubly to 
singly ionized Hehum ratio (He^+ZHe"*"). The computed emissivity coefficient of the nebular 
gas was scaled as an initial approximation to the total flux at the Earth by deriving the 
emitting volume absolute H/3 flux F(H/3), dereddened using chp — 1.475 Eb-v- If the nebula 
was not entirely contained within the aperture used for the observations, the continuum was 
scaled by an appropriate geometrical factor {4A/tt6'^ or 4A/ 710^93 where A is the area of the 
aperture in square arcscconds). Because the value of F(H/3) is very uncertain the continuum 
model was then re-scaled to match the observed flux. 

For NCC 2371, Abell 78, and NGC 1535, the estimated nebular continuum fluxes are 
~ 1%, < 0.1%, and < 0.01%, respectively, of the observed flux between 1400-1500 A. 
These fluxes do not signiflcantly affect the modeling of the stellar spectra. However, the 
contribution of the nebular continuum emission to the UV spectra of IC 4776 is signiflcant, 
and is responsible for the shape of the observed spectrum at the longer UV wavelengths. 
This can be seen in Fig. 7, which shows the UV spectrum of this object along with our 
stellar and nebular continuum models. 



-16- 



4.4. Reddening 

In Table 6, we list reddenings determined from the FUV/UV spectral slope, from our 
measured H I column densities, and from literature. We discuss each in turn. 

For even the coolest objects of our sample, wavelengths above 912 A lie in the Rayleigh- 
Jeans tail of the spectral energy distribution. This makes the slope of the FUV/UV contin- 
uum insensitive to Tg// and mainly dependent on the reddening toward the object. We have 
thus used the slope of the FUV/UV spectra to constrain E^-y toward three of our objects, 
for which the nebular continuum contribution is negligible. In doing so, we have used the 
reddening law of Cardelli et al. (1989) assuming Rv — 3.1, the standard Galactic value. 
Although their work was originally valid only for wavelengths longer than 1250 A, ongoing 
work has shown that the law may be safely extended through the FUSE range (G. Clayton 
2001, private communication). In the case of IC 4776, the nebular continuum contamination 
in the UV (§ 4.3) hinders this method. However the FUV spectrum is primarily stellar, and 
the determination of -Eb-v was based on this range. The sum of the stellar and nebular 
components, and the derived -Eb-v, fit well with the overall wavelength range (900-3000 A). 

Our values determined from the UV slope agree well with those determined from liter- 
ature values of the logarithmic extinction (at 11/5) using the relation = 1.475 £'b-v. In 
all cases, the reddening is very small, and does not affect the results of the stellar modeling. 

Finally, we list the reddenings implied by our measured column densities of H I (§ 4.2) 
using the relationship (A^(H I)/ Eb-y) = 4.8x10^^ cm~^ mag~^ (Bohlin et al. 1978), which 
represents typical conditions in the ISM. This results in significantly higher reddenings than 
the other two methods, indicating that the derived column density includes a significant 
amount of circumstellar H I, which apparently has a smaller dust-to-H I ratio than that of 
the ISM. This is similar to our findings for other CSPN (Herald & Bianchi 2002, 2004). 

5. RESULTS and DISCUSSION 

As discussed in § 4.1, values of log (7 were taken or estimated from previous works based 
on optical spectra (photospheric lines). From our spectral fitting process we determined the 
distance-independent parameters Tg//, Rt, and Voo- Once a distance is adopted, was 
derived by scaling the model flux to the observations, and then L and M were determined. 
Distances to our objects found in the literature derived by various methods show large 
spreads {e.g., for NGC 1535, 0.8 — 3.1 kpc — Sabbadin et al. 1984). Most are determined 
from nebular relations, such as that of the nebular radius to ionized mass. CiarduUo & 
Jacoby (1999) have shown that, in comparison with distance determinations made from 
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more reliable methods such as spectroscopic parallax, statistical methods typically are off by 
a factor of 2 or more. With this in mind, our adopted distances from literature are shown 
in Table 7, along with the resulting parameters. 

The results of comparing the effective temperatures and luminosities of our sample to 
the evolutionary tracks of Vassiliadis & Wood (1994) are shown in Table 8. They include the 
core mass (M^, essentially the current mass), the initial mass of the progenitor (Mj„j(), the 
evolutionary age {rg^oi) and the gravity derived using Mc and the radii (Table 7). NGC 1535, 
the H-rich object, was compared with the H-burning tracks, and lies along the constant- 
luminosity phase. We used the He-burning tracks for NGC 2371 and Abell 78, which both 
fall on the bend between the constant-luminosity phase and the WD cooling sequence, as 
appropriate for their transitional nature. At the temperature of IC 4776, the H-burning 
tracks extend down to luminosities of logL > 3.5 Lq (corresponding to the M^nit — 1-0 M0 
track). The upper hmit of our luminosity of IC 4776 (logL = 3.21q;^5 Lq) does not overlap 
with the H-burning tracks. The parameters of IC 4776 do, however, fall on the He-burning 
tracks, which we have used in this case. Which tracks were used are listed in the Table. 

In most cases, masses calculated using our derived radii and initially adopted gravities 
( log ^ = 5.5 for Abell 78 and NGC 2371, log g = 4.3 for NGC 1535 and IC 4776 — § 4.1.2) are 

unreasonably low {i.e., < 0.5 Mq). Gravities computed using the track masses are higher 
than our adopted values in all cases. The gravities of NGC 2371 and Abell 78 (log (7 = 
6.3 and 5.7, respectively) are good agreement with PC 1159 stars of similar temperatures 
from the 16-objcct sample of Miksa ct al. (2002). Likewise, the gravities of NGC 1535 
and IC 4776 {logg = 4.9 and 5.1) are reasonable for CSPN of these temperatures {e.g., 
Mendez et al. 1988). Models with these higher gravities resulted in only slight changes to 
the stellar parameters, and within the errors quotes in Table 7. Therefore we conclude these 
higher gravities are more appropriate for our sample than those originally adopted from the 
literature. 

For most of our objects, the dynamic ages (Table 5) are 2-4 times lower than the 
evolutionary ages. This is not unexpected, as the dynamic age is a lower limit to the post- 
AGB age, because the nebular expansion increases during the early post-AGB phase and then 
level off as the nucleus fades (see Sabbadin 1984, Bianchi 1992 and references therein). The 
exception is Abell 78, which has a larger dynamic age than its evolutionary age {rdyn — 15 kyr 

vs. Tevol = 10). 

Vink et al. (2001) presented wind models that yield mass-loss rates for given stellar 
parameters, taking into account opacity shifts in the wind due to different ionization struc- 
tures at different temperatures. They assume "normal" abundances - i.e., not He or CNO 
enriched. Although their prescriptions were derived from studies of the winds of (massive) 
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O and B stars, we applied their prescription using the derived parameters {Teff, L, Voo, and 
log (7) of our low-mass objects (Table 7), using solar metallicity. The results arc listed in 
Tabic 7. The predicted mass-loss rate of NGC 1535 (which we fit with normal abundances) 
is the same as our measured value. For NGC 2371 and Abell 78 (the H-deficient objects), 
the predicted mass-loss rates are many magnitudes lower than our values. This is probably 
due to the chemically enriched winds of these objects having much higher opacities. For 
IC 4776, a predicted mass-loss rate is a magnitude lower than our measured value results. 
This may be an indication that the wind of IC 4776 has some chemical enrichment, another 
indication it may be an He-burner. 

The mass-loss rates of the H-dcficient objects are about a magnitude larger than the 
H-rich object, and have higher terminal velocities. The more powerful winds of the hotter 
objects are a consequence of the hydrogen deficiency of the atmosphere, which increases the 
opacity and thus the efficiency in converting radiative momentum flux into wind momentum 
flux. A measure of the efficiency of the wind at converting radiative momentum flux to 
wind momentum flux is the "performance number" r) — VcoMc/L (Springmann 1994). The 
performance number basically measures how many times, on average, each photon scatters 
in the wind. Performance numbers above unity (the "single scattering limit") are typical 
of stars with chemically enriched winds such as WR objects. The performance numbers for 
our sample are also fisted in Table 7, and show the H-rich object having the lowest 77, while 
the two H-deflcient objects are above the single-scattering limit. IC 4776 has a performance 
number of unity, again indicating its wind may be chemically enriched. 



6. CONCLUSIONS 

We have performed a FUV/UV spectral analyses on four CSPN which display com- 
bined absorption and wind line spectra, a sign that the winds of these objects are in the 
process of switching off. Two objects are H-deficient (NGC 2371 and Abell 78), one is H-rich 
(NGC 1535), and for one our analysis is inconclusive (IC 4776). Thus both channels of the 
CSPN phase are represented. We derived the stellar wind parameters of the sample PN by 
modeling their FUV and UV spectra, where strong wind signatures are seen. NGC 1535 and 
Abell 78 have been previously analyzed in the optical using plane-parallel models(Werncr 
et al. 2003, Mendez et al. 1988), and the latter also with wind models (Koesterke & Werner 
1998). Our parameters are in good agreement with the previous results. 

The most striking feature in all the FUV spectra is the O VI doublet. In massive 0-type 
stars, this ion originates from X-rays due to shocks in the radiatively driven winds (Bianchi 
& Garcia 2002). We are able to fit this feature in our sample CSPN without the inclusion 



-19- 



of X-rays due to their higher effective temperatures. 

Two objects have T^ff ~ 65 kK, with M ~ 10"^ yr^^ NGC 1535 has been classified 
as a Hydrogen- rich 0(H) star based on its optical spectrum, and we are able to match its 
FUV/UV spectrum using solar abundances. Prom the spectrum of IC 4776, we cannot 
constrain the abundances, but comparison with evolutionary tracks suggest this object is 
a He-burner. Our models for these objects show a jump in the ionization structure of the 
CNO elements (and S and P) between 55-70 kK. It is in this temperature range that Tinkler 
& Lamers (2002) find that H-rich CSPN exhibit a large jump in M by an factor of 10-100, 
which they suspect is due to such a bi-stability jump. This supports the idea of Vink ct al. 
(2001) that a bi-stability jump at such high temperatures may be due to CNO elements. 

The two H-deficient objects Abell 78 and NGC 1271 are significantly hotter (Tg// > 
110 kK), with M ~ lO"'^ M© yr'^ Our derived parameters for the [W01]-PG 1159 star 
Abell 78 generally agree with those of past analyses. We derive similar parameters for 
NGC 2371, and suggest that it is a [WR]-PG 1159 also, but its wind is more ionized and 
shows no O V features (a [WO0]-PG 1159?). We find evidence of iron deficiency in both 
of these objects, supporting the findings of Werner et al. (2003) in Abell 78 and of Miksa 
et al. (2002) in PG 1159 stars. They he on the bend in the theoretical evolutionary tracks 
between the constant luminosity phase and the WD cooling sequence, having post-AGB ages 
of 10-15 kyr. For [WC] stars, it seems that as the star evolves away from the post-AGB phase 
through the [WC] sequence, the temperature and terminal velocity of the wind increase as 
the wind density decreases (Acker & Neiner 2003). 

Except for Abell 78, the post-AGB ages predicted by the evolutionary models arc typi- 
cally 2-4 times lower than their kinematic ages. However, kinematic ages are lower limits to 
the post-AGB age (as the expansion is slower in the initial phase) thus the actual discrepancy 
may be smaller. 

Our FUV/UV analysis has provided wind parameters for H-rich and H-deficient CSPN 
at the stage in post-AGB evolution where the winds are fading. However, because a precise 
determination of the mass is not possible, the objects do not necessarily represent the same 
evolutionary sequence. This work has also provided information on the interstellar and 
circumstellar environment from our measurements of the column densities and temperatures 
of H I and H2 along the sight-lines. Our determinations of the H I column density and 
£'b_v imply that the relationship between these two quantities in the circumstellar (PN) 
environment of CSPN differs from that of the ISM, having lower dust-to- gas ratios (probably 
due to the destruction of dust by the radiation field). The high resolution of the FUSE data 
allow us to detect hot H2 associated with the nebulae. For all four objects in our sample, a 
single component of H2 gas at typical ISM temperatures (i.e., T ~ 80 kK) was not adequate 
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to fit the absorption spectrum. A second, hotter {T ~ 300 kK) component was necessary, 
which we assume to be associated with the circumstellar environment. With the advent of 
many FUSE observations of CSPN, it is becoming apparent that hot circumstellar H2 is not 
uncommon CSPN at quite different evolutionary stages: from old FN with white dwarf nuclei 
(Herald & Bianchi 2002) to young, compact FN (Herald & Bianchi 2004). The objects in this 
paper lie either along the constant luminosity section of the post-AGB evolutionary tracks, 
or on the transition bend to the WD cooling sequence, and thus represent intermediate 
stages to the previously mentioned cases. Given the intense UV radiation fields emitted by 
the CSPN, it is likely that the nebular H2 exists in clumps, shielded by neutral and ionized 
hydrogen, as appears to be the case in the Helix nebula (Speck et al. 2002). 
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Telescope Science Institute (STScI) . STScI is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS5-26555. 



A. APPENDIX: MODEL ATOMS 

Ions and the number of levels and superlevels included in the model calculations are 
hsted in Table 9. The atomic data come from a variety of sources, with the Opacity Project 
(Seaton 1987; Opacity Project Team 1995, 1997), the Iron Project (Pradhan et al. 1996; 
Hummer et al. 1993), Kurucz (1995)^ and the Atomic Spectra Database at NIST Physical 
Laboratory being the principal sources. Much of the Kurucz data were obtained directly 
from CfA (Kurucz 1988, 2002). Individual sources of atomic data include the following: 
Zhang & Pradhan (1997), Bautista & Pradhan (1997), Becker & Butler (1995), Butler et al. 
(1993), Fuhr et al. (1988), Luo & Pradhan (1989), Luo et al. (1989), Mendoza (1983, 1995, 
private communication), Mendoza et al. (1995), Nussbaumer & Storey (1983, 1984), Peach 
et al. (1988), Storey (1988, private communication), TuUy et al. (1990), and Wiese et al. 



^See http: / /cfa- www.harvard.edu/amdata/ ampdata/amdata.shtnil 
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(1966, 1969). Unpublished data taken from the Opacity Project include: Fe VI data (Butler, 
K.), Fe VIII data (Saraph and Storey) and Fe IX and Fe X data (C. Mendoza). 
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Fig. 1. — FUV/UV spectra of our sample. Data are from different instruments (with different 
resolutions - see Table 2). The spectrum of IC 4776 appears flatter than the others because 
of the contribution from the nebular continuum (see § 4.3). Prominent stellar features 
are marked with dark/blue labels, and interstellar/nebular features by light/red features. 
Regions of bad data are blanked out. 
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Fig. 2. — The ionization structure of iron for models spanning the temperature ranges of 
our objects, as a function of distance from the stellar surface. The upper panels are models 
with temperatures of our cooler objects {Tgff ~ 60 — 70 kK), where Fe VI-Fe VII dominate 
throughout the wind. The lower panels models in the temperature range of our hotter 
objects {Teff — 100 — 130 kK), where the higher ionization stages arc dominant. These 
figures illustrate the sensitivity of the ionization structure of iron to T^ff. 
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Fig. 3. — NGC 1535: Iron spectrum as a diagnostic. Observations of NGC 1535 are shown 
(black) along with a model with T^ff ^ 75 kK, M ~ 3x10^® Mq yr^^ (green/light gray). 
This model displays an Fe VII spectrum stronger than the observations, indicating a lower 
mass- loss rate and/or a different temperature is required. Our final model, with T^ff ^ 
65 kK, M ~ lxlO~^ Mq jt~^ (red/heavy gray) provides a better fit. 
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Fig. 4. — Model fits. The observations (black) have been unreddened using our values for 
Eb-y derived from the UV slope, and our models (green/thick gray) are also shown. In the 
case of IC 4776, the model nebular continuum (blue/dark dashed), and the stellar model 
(pink/light gray) are summed together (green/thick gray). The prominent wind features 
are marked with black labels, and strong nebular and interstellar features are marked with 
gray/red labels. The spectra have been convolved with a 0.75 A Gaussian for clarity. 
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Fig. 5. — Fitting H2: The FUSE spectrum of NGC 1535 is shown (black) along with a 
stellar model with successive hydrogen absorption components applied (Table 4). The dashed 
red/dark gray line shows the stellar model with absorption from the "cool" (T = 80 kK) 
H I and H2 components applied, and the green/light gray line with absorption effects of the 
"hot" H2 component (T = 400) H2 applied to that. Features attributable to the hot and 
cool components, as well as strong interstellar absorption features, are labeled. The cool 
hydrogen absorption model is not sufficient to fit all the features, suggesting a hot H2 gas 
associated with the nebula. 
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Fig. 6. — Fitting Lya of NGC 1535: The Lya profile in the STIS data (black) and our stellar 
model with the absorption effects of different amounts of H I applied (T = 80 kK): column 
densities of logA^ = 20.8, (green/light gray) and 21.0 cm~^ (dashed red/dark gray). The 
models bracket the observations, so the sight-line column density can be well constrained to 
logA^ = 20.8 ±0.1 cm-2. 




Fig. 7. — IC 4776: The observations are shown (gray) along with our stellar (dotted) and 
nebular continuum models (dashed) . The sum of both models is also shown (solid) , demon- 
strating the flux longwards of Lya (1215 A) of this object is signiflcantly affected by the 
nebular continuum. We have applied our molecular hydrogen absorption model (§ 4.2) to 
the stellar model, and our atomic hydrogen absorption model to both the stellar and nebular 
models. In the FUV, the observations and models have been convolved with a 2 A Gaussian 
for clarity. 
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Fig. 8.— O VI AA1031.93,1037.62 for NGC 2371, in velocity space around 1031.93A. The 

spectra has been blueshifted to account for the measured radial velocity of the star of 
Vr — +20.6 km s^^ (Table 1). Marked is a narrow emission line at ~ 20 kms^^, probably 
corresponding to the blue component of the doublet. We have also marked the corresponding 
position of the red component, assuming it is velocity shifted by the same amount. The red 
component seems to be present, but the absorption from H2 prevents positive identification. 
Kaler et al. (1993) noticed similar narrow emission lines on top of the broad P-Cygni profiles 
in optical O VI features. 
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Table 1: The Program Objects 

Name PNG R.A. Dec 'i^ PN Diameter 

(J2000) (J2000) (kins--^) (") 



NGC2371 
Abell 78 

NGC 1535 
IC 4776 



189.1+19.8 

081.2-14.9 
206.4-40.5 
002.0-13.4 



07 25 34.71 
21 35 29.8 
04 14 15.9 
18 45 50.57 



+29 29 26.20 

+31 41 40 
-12 44 21 
-33 20 32.94 



+20.6 + 2.7 
+17+10.0 
-3.2 + 1.4 
+18.9 + 0.7 



48.9x30.6*^ 

107'^ 
33.3x32.1^ 
8.5x4.0'' 



a: Cahn & Kaler (1971) 

b: Tylenda et al. (2003), two dimensions given. 



Table 2: Utilized Spectra 



Star 


Instrument 


Data- 


Date 


Resolution 


Aperture 


Range 


Scale 






set 






(") 


(A) 


Factor 


NGC 2371 


FUSE 


P1330301 


02/26/00 


~ 0.05 A 


30x30 


915-1180 






lUE 


SWP04883 


04/07/79 


5-6 A 


10x20 


1180-1975 


xl.3 




lUE 


LWR04210 


04/07/79 


5-6 A 


10x20 


1975-3345 


xl.3 


Abell 78 


BEFS 


BEFS2190 


12/01/96 


~ 0.33 A 


20 


915-1222 






lUE 


SWP19906 


05/05/83 


~ 0.2 A 


10x20 


1122-1975 


xl.2 


NGC 1535 


FUSE 


P1150808 


10/05/01 


~ 0.05 A 


30x30 


915-1180 






STIS+E140M 


O64D04010 


03/01/01 


~ 0.05 A 


0.2x0.2 


1180-1700 




IC 4776 


FUSE 


P1330501 


05/21/00 


~ 0.05 A 


30x30 


915-1180 






lUE 


SWP16504 


03/11/82 


5-6 A 


10x20 


1180-1975 


xl.25 




lUE 


LWR12764 


03/11/82 


5-6 A 


10x20 


1975-3345 


xl.25 
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Table 3: Model abundances ( Xq = solar abundance) 



Parameter 


Xh 


Xhc 


Xc 


Xn 


Xo 


Xsi 


^5 


Xpe 


NGC 2371 




0.54 


0.37 


<o.oi? 


0.08 


Xo 


Xo 


Xo? 


Abell 78 




0.54 


0.37 


0.01 


0.08 


Xo 


Xo 


<0.03Xo 


NGC 1535 


Xo 


Xo 


Xo 


Xo 


Xo 


Xo 


Xo 


Xo 


IC 4776 




Xo 


Xo 


Xo 


Xo 


Xo 


Xo 


Xo 



?: Value uncertain, see text. 



Tabic 4: H2 and H I parameters 



Star 


Component 


logiV 


T 


Note 


NGC 2371 


H I(IS+circ) 


21.0 


80 


Not fit 




H2 (IS) 




80 






H2 (circ) 


17.0+°-? 


300 




Abell 78 


H I(IS+circ) 


21.1^^2 


80 


Prom Lya 




H2 (IS) 


i9.7t^:? 


80 






H2 (circ) 


16 4"°-'^ 


300 




NGC 1535 


H I(IS+circ) 


20.8±0.1 


80 


Prom Lya 




H2 (IS) 




80 






H2 (circ) 


16.4l°i 


400 




IC 4776 


H I(IS+circ) 


21.4l°i 


80 


Prom Ly/5 




H2 (IS) 


1 5 7"°-'^ 
-'-'J- '+0.3 


80 






H2 (circ) 




300 





Table 5. Nebular Parameters and Estimated Kinematic Ages 



Name 


rie 
(cm-3) 


Te 

(K) 


(ergs cm~^ s~^) 


He/H 


He2+/He+ 


Ve.p{0 III) 

(kms"^) 


^adopt 

(") 


(pc) 


'^dyn 

(kyr) 


NGC 2371 


1260" 


14,100'' 


-10.99'' 


o.ir 


o.ess** 


42.5'' 


40 


0.15 


3.3 


Abell 78 


790" 


20,800'' 


-12.04'' 


0.110'' 


1.61'' 


27/ 


107 


0.42 


15 


NGC 1535 


21009 


10,8009 


-10.45'' 


0.096'' 


0.126'' 


20/ 


33 


0.13 


6.3 


IC 4776 


12,000'' 


8600''''* 


-10.73'' (-10.43^ 


0.085'' 


0.005'' 


22^' 


6 


0.06 


2.5 



References. — (a): Feibclman (1997) (b): Cahn ct al. (1992) (c); Kalcr ct al. (1993) (d): Sabbadin (1984) (c): 
Medina & Pena (2000) (inner knots, upper limit) (f): Weinberger (1989) (g): Tylenda et al. (1991) (h): Phillips 
(1998) (i): adopted, see text, (j): Gesicki & Zijlstra (2000) 

I 

CO 

O) 

I 



-37- 



Table 6: Reddening parameters 



Star -E-B-v E-q-y -E-b-v 

(UV slope) (H I) 

(mag) (mag) (mag) 



NGC 2371 


0.070±0.02 


0.07" 




Abell 78 


0.100±0.02 


0.12^ 


0.21 


NGC 1535 


0.055±0.02 


0.07^ 


0.13 ±0.03 


IC 4776 


0.090±0.02 


0.07^ 


0.52l°i? 



a: Kalcr ct al. (1993) 
b: Cahn et al. (1992) 



Table 7: Derived Stellar Parameters and Adopted Distances 



Star 


Teff^ 






D 


Rj 


log it 


log 


log 






(kK) 


(Ro) 


( km s"'^) 


(kpc) 


(Ro) 


(Lo) 


(Mq yr-1) 


(Mq yr-i) 




NGC 2371 


135tJ^ 


15t^° 


3700 ± 200 


1.5" 


0.09 


^•^"3-0. 20 


-7.11 ±0.30 


-10.6 


5.0 


Abell 78 


113 ±8 


07+20 


3200 ±50 


1.6" 


0.19 


q 70 + 0.10 

' '-'-0.13 


7 qq+0.36 
1 .00_Q 


-9.3 


1.4 


NGC 1535 


66 ±5 


200+^f 


1950 ±50 


i.e*" 


0.43 


3.60 ±0.13 


-8.10 ±0.30 


-8.1 


0.2 


IC 4776 


601^° 




2300 ± 200 


3.9" 


0.35 


o.^U_o.i5 


-7.85 ±0.30 


-8.8 


1.0 



(a) : Cahn et al. (1992) 

(b) : Sabbadin et al. (1984) 



f: Prom our spectral modeling. 

p. Predicted mass-loss rate for our derived parameters T^ff, L, logg, Voo, using prescription of Vink et al. (2001). 



Table 8: Derived Stellar Parameters From Evolutionary Tracks of Vassiliadis & Wood (1994) 

Star Track Mc M^mt logs 

(Me) (Me) (kyr) (cms-^) 

NGC 2371 He 0.63 2 ~ 13 6.3 

Abell 78 He > 0.63 > 2.0 < 10 5.7 

NGC 1535 H 0.58 1.2 ~ 15 4.9 

IC 4776 He 0.57 1.0 ~8 5.1 
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Table 9: Levels and super levels for model ions 

Element I H HI fv V VI VH VHI IX X XI 

~H. 20/30 TJl 

He 40/45 22/30 1/1 



c 


30/54 13/18 


1/1 






N 


29/53 


13/21 


1/1 







29/48 


41/78 


13/19 


1/1 


Si 


22/33 


1/1 






P 


36/178 


16/62 


1/1 




S 


51/142 


31/98 


28/58 


1/1 



Fe 51/294 47/191 44/433 41/254 53/324 52/490 43/210 1/1 



